Universal 1+2-body bound states in planar atomic wave guides 
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Shallow heteronuclear trimers are predicted for mixtures of two atomic species strongly trapped in 
a quasi two-dimensional atomic wave guide. The binding energies are functions of the 2D-scattering 
length and of the mass ratio and can be thus tuned by various ways. These universal trimers 
are composed of two identical non interacting particles and of a third particle of the other species. 
Depending on the statistics of the two identical particles, the trimers have an odd (fermions) or even 
(bosons) internal angular momentum. These results permit to draw conclusions on the stability issue 
for the quasi-2D gaseous phase of heteronuclear dimers. 
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INTRODUCTION 



Actual experiments in ultracold physics open the possi- 
bility of trapping and coolingatomic mixtures composed 
of different atomic species in quasi two-dimensional 
(quasi-2D) geometries Q. These dramatic progress 
open exciting possibilities for exploring experimentally 
few- and many-body properties in these systems. For 
example, the Bcrczinski-Kostcrlitz-Thouless transition in 
bosonic 2D gases is under active experimental studies @- 
p. The occurence of Confinement Induced Resonances 
[10l4l2| i.e. for quasi-2D wave guide, resonances in two- 
body scattering induced by the transverse confinement, is 
another peculiar property in these systems which repre- 
sents a way to achieve strongly correlated regimes. There 
are also promising possibilities using fermionic quasi-2D 
gas to reach exotic phases 13}. More recently, a new 
crystalline order has been predicted to occur for two- 
species quasi-2D fermionic gases with a sufficiently large 
mass ratio |14] . This last reference suggests also the idea 
of tuning the mass ratio by imposing an optical lattice 
acting on one species only. This opens fascinating possi- 
bilities for achieving systems with very high mass ratio. 

Motivated by these recent progress, we consider in this 
paper mixtures of two atomic species 'strongly confined' 
along the z-direction by a planar atomic wave guide of 
frequency uj z , without any external confinement along the 
x and y directions. Here, the expression 'strongly con- 
fined' means that for large interparticle separations {i.e. 
for distances which are very much larger than the trans- 
verse harmonic oscillator length denoted by a z ), atoms 
are frozen in the ground state of the ID harmonic os- 
cillator along z. In this situation, the system can be 
considered as quasi-2D and long range properties (for in- 
terparticle distances very much larger than a z ) acquire 
the peculiar aspects typical of a 2D-world. The results of 
the paper can be summarized as follows: (i) we extend 



the formulation of the 2D zero-range potential (ZRP) ap- 
proach which is usually used in the configurational space 
[15l420| | to the momentum representation which is more 
natural for a translationally invariant problem. This way 
the equivalence between the ZRP model and the diagram- 
matic [2l| or the i-matrix approaches [12] is much more 
transparent, (ii) We predict the existence of '1 + 2-body' 
universal shallow bound states composed of two identical 
atoms of mass M and of a third atom: an 'impurity' of 
mass rrii. There is no interaction between the two iden- 
tical atoms of mass M while each one interacts (in the s- 
wave channel only) with the impurity. For each statistics, 
the bound states are labeled by a radial quantum num- 
ber and an internal angular momentum. This later one 
is even when identical atoms are bosons and odd when 
they are fermions. The binding energies are function of 
the two-body scattering length, of the oscillator length 
and of the mass ratio, thus it can be tuned by various 
ways. If the two identical atoms obey the Fermi statis- 
tics, we find that for increasing values of the mass ratio 
x = M/rrii, the first trimer appears at x — x\ ~ 3.33 for 
a p-wave internal angular momentum, (iii) Neglecting 
inelastic processes toward deep molecular bound states 
with spatial extensions of the order of the range of inter- 
atomic forces (which is very much smaller than a z ), the 
low energy two-body collisions between 2D shallow het- 
eronuclear dimers composed of fermions are elastic for a 
mass ratio x < X2 — 18.3. (iv) In the case where identical 
atoms are bosons, at least one trimer exists for all values 
of the mass ratio. Moreover, binary collisions for dimers 
composed at least of one boson can always populate a 
shallow trimer state. 

The structure of the paper is as follows. First, we recall 
the low energy two-body scattering properties for atoms 
strongly confined in a 2D plane by a harmonic trap along 
the third dimension. We then introduce the ZRP model 
in the momentum representation, where the contact con- 
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dition for a vanishing distance between interacting par- 
ticles is replaced by an integral equation. In the second 
part, we explain how this formulation can be extended 
to few-body problems. In the third part, we show that 
depending on the various parameters (statistics, mass ra- 
tio, 2D-scattering length) the three-body wave equation 
can support one or more bound-states. We then discuss 
about the implications of these findings. In particular we 
focus on the stability issue for the quasi-2D heteronuclear 
dimer gas. 



I. MOMENTUM REPRESENTATION OF THE 
ZERO-RANGE POTENTIAL MODEL 

In this first part, we introduce the ZRP model in the 
momentum space. We then consider a system composed 
of two interacting particles of respective masses M and 
m% (with the reduced mass /i) interacting in the s-wave 
channel only. The particles are trapped by a harmonic 
potential along the z-axis with the same trap frequency 
(u z ) for both species. We then define the harmonic os- 
cillator length as: 




(1) 



According to the 'strong 2D confinement' hypothesis, 
we suppose that for large interparticle separations (very 
much larger than a z ), the wavefunction factorizes in the 
product of a wavefunction of the 2D coordinates (2D- 
wave function) and of two ground state eigenfunctions of 
the ID-harmonic trap for each atomic z-coordinate. In 
what follows, a z represents a short range scale (in the 2D- 
world, this length is formally equivalent to the range of 
the true interatomic forces in 3D) and we consider only 
the quasi-2D behavior: atomic momenta and positions 
are in the 2D xy plane. The two-body scattering states 
for an incident plane wave of wave vector ko are given in 
the low energy regime (koa z -C 1) by: 



(k|f ko ) = (2n) 2 S(k - k 
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where fw(ko) is the 2D scattering amplitude [23 
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In Eq. ([3]), 7 = 0.577 • • ■ is the Euler's constant and <Z2d 
is the so-called 2D-scattering length. For a pairwise in- 
teraction in 3D characterized by a scattering length a 3 rj , 
d2D is a function of the characteristic length of the trap 
a z and also of the three dimensional scattering length 
«3D Q: 



a 2 D = 



VeTexp(J) 



exp 



a z y/lT 



2 a 



3D 



(4) 



and the constant J in Eq. @ is [H, [24 



J = 



du 



= -1.36 •• 



(5) 

Eqs.([2l [3]) can be deduced by using a 2D ZRP model, 
which by construction is valid for large interparticle spac- 
ing with respect to a z 1 1 81-4201] . Following the spirit of the 
Bethe-Peierls approach[25f7we recall briefly the formal- 
ism in the configurational space for a stationary state | \?) 
of energy E in the center of mass frame. The pairwise 
interaction in the non-interacting Schrodinger equation 
is replaced by a S(r) source term characterized by an 
amplitude denoted in what follows by K: 



(r|tf) = —5(r)K. (6) 



Solutions of Eq. ([5]) are linear combination of a singular 
and of a regular function. For a regular function sat- 
isfying the boundary conditions associated with a given 
situation (for example an incident plane wave in the case 
of a scattering problem), the interacting wavefunction is 
obtained by determining the balance between these two 
functions i.e. by choosing a specific value for the source 
amplitude K. In the present case and analogously to 
what happens for point-charge in 2D-electrostatic, the 
source term leads to a logarithmic singularity for vanish- 
ing interparticle separations r — > 0. The correct balance 
between the singular and the regular contributions is de- 
duced by imposin g th e following behavior on the wave- 
function as r jl5ll2C|: 



(r|*) 



-Kin 



\«2D 



O(r) 



(7) 



This contact condition allows one to recover the scatter- 
ing states in Eq. ([2]) exactly. It is worth pointing out 
that the logarithmic singularity is purely formal because 
for r of the order or smaller than a z , the 2D zero range 
model ceases to be valid and 3D physics becomes impor- 
tant. However, this is the key property which permits 
to implement the effect of the interaction for interparti- 
cle separations very much larger than a z . The systems 
which are studied in this paper are translationally invari- 
ant along the xy plane. Hence, solving the stationnary 
Schrodinger equation is a much more simple task in the 
k-representation than in the configurational space. Con- 
sequently, we formulate below the 2D ZRP model in the 
momentum representation. For this purpose, we closely 
follow the method already introduced in Ref. [HJ for ar- 
bitrary resonant Z-wave interaction in 3D. First, we in- 
troduce a source term with a delta of finite range e: 



(k|<5 £ ) = exp 



(8) 



which is such that: lim e _j. ( r l^e) — S(r). Keeping the 
same notations as in Eq. ©, the stationary states \^> e ) 
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at energy E verify: 

'h 2 k 2 



^ £ 1 <k|* £ ) = — <k|«5 e )ft. (9) 



and the stationary states in the zero-range limit are: 



I*) = lim |* e 

£->0 



(10) 



Let us consider the 2D Green's function of Eq. ([6]) (solu- 
tion corresponding to 1Z — 1) taken at energy — fi, 2 A 2 //i < 
where A is any real and positive number. This func- 
tion is Ko(Ar) and behaves as — ln(Are 7 /2) + 0(r) for 
r — > 0. Using this property, the contact condition in 
Eq. (J7J can be written in terms of a subtraction between 
the eigenstate of the ZRP model and this Green's func- 
tion as r — > 24]. In the k- representation the contact 
condition can thus be written as: 



lim 



d 2 k 

(2n) 2 



(kltf A - 



2ir(k\5 e )1l 
k 2 + A 2 



K 
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where the scattering amplitude at negative energy is: 

1 



/a D (iA) 
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A > 0. (12) 



We emphasize that the interacting wave function ob- 
tained by using Eq. (fTTjl do not depend on the parameter 
A: this is the so-called A-freedom introduced previously 
through the A-potential in Ref. [18j. As an example, we 
illustrate this method in the case of a two-body collision 
for an incident plane wave of momentum ko in the center 
of mass frame. For E > 0, using the usual prescription 
+ i0 + corresponding to an outgoing cylindrical 
wave, Eq. @ can be transformed into: 



(k|* e ) = (2^) 2 <5(k-k ) + 



irh 2 



h?k 2 

2/j 



E -i0+ 



(13) 



Then, inserting Eq. (IT31 in the contact condition 
Eq. (fTT|) . one obtains 1Z = /2d(^o) given by Eq. ([3]) 
VA > 0. Therefore as expected, performing the e — >• 
limit in Eq. (|13[) one recovers the 2D s-wave scattering 
states of Eq. ([2]). The 2D ZRP model supports also the 
existence of a two-body bound state characterized by a 
binding energy Ed im = —E where: 



2(i 



with g<jir 



e 7 a 2D ' 



(14) 



and a wavefunction \4>b) given in the r- representation by: 

9dim 



(r|0i 



-K {qdimr) 



(15) 



This bound state has a spatial extension of the order of 
a 2 D and can thus be considered as shallow in the limit: 



»2D > a z 



(16) 



Eq. (|T6|) defines the quasi-2D resonant regime where the 
extension of the 2D dimer is large with respect to a z and 
is analogous to the 3D resonant regime where the 3D scat- 
tering length is large in absolute value with respect to the 
typical radius of interatomic forces. The length a z gives 
the order of magnitude of the radius of the 2D effective 
interaction [TTj so that the quasi- 2D resonant regime de- 
fines also the regime where the zero-range approximation 
is well justified. For instance, the scattering amplitude 
in Eq. ^ has a maximum modulus for k = qa\m and in 
the quasi-2D resonant regime, the maximum occurs for a 
2D momentum which is very small with respect to 1 j a z 
i.e. for a collisional energy such that the 2D approxi- 
mation remains valid. For a fixed value of the collisional 
momentum, the maximum can be reached if one tunes 
the dimer binding wavenumber gdim (and thus the 2D 
scattering length) by modifying the 3D-scattering length 
and/or the trap frequency oj z [see Eq. J3J)]: this is the 
so-called 2D Confinement Induced Resonance first found 
in Ref. (l0|. As a consequence of Eq. (HJ), the condition 
in Eq. (fTB)) corresponds to the regime of negative 3D- 
scattering lengths where also |o3d| is much smaller than 
a z . In this paper, we suppose that Eq. (|16p is satisfied 
for heteronuclear atomic pairs. 

Note that in the case where an optical lattice in the 
xy-plane acts for example on the heavy atoms of mass 
M in the small filling limit as proposed in Ref. [14||, we 
assume that <z 2 d is much larger than the lattice period. 
In this limit the mass M can be replaced by the effective 
mass in the 2D-Schr6dinger equation [261 ] . 



II. FORMULATION OF THE 1 + 2-BODY 
PROBLEM 

In this second part, we show how the ZRP model can 
be implemented in few-body systems and in the momen- 
tum representation by considering a particular '1 + 2- 
body' problem. Our system is composed of two identical 
atoms of mass M each interacting only with a single atom 
impurity of mass rrii via the ZRP model. Hence in or- 
der to fulfill the strong 2D confinement condition's, we 
require implicitly that the collisional energies of the prob- 
lem are very much smaller than the characteristic energy 
of the transverse trap (haj z ). 

Furthermore, we assume that the quasi-2D resonant 
regime defined in Eq. (|16[) is always achieved for het- 
eronuclear atomic pairs. We introduce the statistical 
parameter 77 which is equal to —1 if the two identical 
atoms are fermions, and to +1 if they are bosons. In the 
fermionic case, the two identical atoms do not interact in 
the s-wave channel as a consequence of the Pauli princi- 
ple. In the bosonic case, we also consider that the identi- 
cal atoms do not interact in the sense that their pairwisc 
interaction is negligible with respect to the interaction 
with the other atomic species. This hypothesis is justi- 
fied in the case where the 3D scattering length between 
the two identical bosons is positive and small in compari- 
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son with a z which implies that their 2D scattering length 
is exponentially small in comparison to a z [see Eq. Q]. 
Consequently, in both cases we take into account the pair- 
wise interaction only between particles of different mass. 
In the ZRP model, for each interacting pair the in- 

teraction is implemented by inserting in the Schrodinger 
equation a source term with an amplitude denoted by 
1Z l± ^ 3 . Consequently, we assume implicitly that configu- 
rations where three particles are in a hyper-radius smaller 
or of the order of a z can be neglected. This hypothesis is 
especially justified in the quasi-2D resonant regime where 
the characteristic extension of the trimers (and thus of 
the hyper-radius) is of the order of a-iD which is very 
large with respect to a z . In the following, we associate 
the three momenta |k;} = {ki,k2,k3} with the particle 
configuration: (1 : M ; 2 : M ; 3 : rtii), furthermore we 
solve the problem in the center of mass frame so that: 



ki + k 2 + k 3 = 0. 



(17) 



For each pair the relative coordinates denoted by 

77; (where i,j,l are in cyclic order) verify: 



k 2 - xk 3 
1 + x 



f)2 



zk 3 - ki 
1 + x 



ki 



(18) 



where x is the mass ratio x — M/rrii. Using these nota- 
tions, the Schrodinger equation for an eigenstate |\& e ) of 
energy E reads in the momentum representation: 



h 2 (k 2 



+ 



h 2 k 2 



2M 2mi 
K^ 3 {r, 2 \6 e )+K 2 ^ 



({k<}|*« 

*7ll<5e>l. 



irh 2 



(19) 



In Eq. (fT9j) . 1Z %± ^ 3 is a function of the total momentum 
of the pair (i,j) and of the momentum of the third par- 
ticle (I). Using the fact that we solve the problem in 
the center of mass frame, TZ l± ^ 3 can be considered as a 
function of the momentum k;. Moreover, the statistics of 
the two identical particles imposes a symmetry between 
the two source terms: 7?. 1± ^ 3 (k) = r/TZ 2 ^ 3 (k) . In order 
to simplify the subsequent equations, we will denote the 
amplitude 7Z 2± ^ 3 (k) by 7Z(k). We also consider in what 
follows situations of negative total energy E only, we thus 
introduce the momentum q: 



2„2 



E 



K 2 q 
2/j 



where q > 0. 



(20) 



The equation satisfied by the function 7Z(k) is deduced 
from the contact condition applied to an interacting 
pair. For example, application of Eq. (jTTJ) for the pair 
(2 : M ; 3 : m$ ) of relative momentum r) 1 gives: 



(2tt)2 



<{ki}|*> - 



27rft(ki 



r\\ + A 2 



ft(ki) 

/2D(iA)" 



For E < 0, Eq. (|2"Tj) can be simplified as: 



^(k) 

/2D(Kfc) 



d 2 U 



(2tt) 2 u 2 



2-k1Z(vl) 



2yk.u + q 2 ' 



(21) 



(22) 



where the variable y is a function of the mass ratio: 
M x 



y 



(23) 



M + rrii 1 + x' 
and also we have introduced the collisional momentum 



W 1 



(24) 



The physical interpretation of Eq. (|24|) proceeds as fol- 
lows. The scattering process between the two interacting 
particles (2, 3) occurs in their center of mass frame at a 
kinetic energy E j which verifies: 



E = E. 



(2,3) 
col 



fi 2 (k 2 



2(M + rrii) 



h 2 

— — k 2 



(25) 



In Eq. (|25p we have used the fact that, as mentioned pre- 
viously, the third atom (1 in this case) does not interact 
with an atom of the pair (2, 3) during the scattering pro- 
cess. Hence, using Eqs. (|24l25p the collisional energy of 
the pair can be written as: 



E. 



2(i. 



<0, 



(26) 



and Kki is the relative momentum of the pair (2, 3) of 
total momentum (— ki). Eq. (|22[) is the 2D analog for 
the 1 + 2-body problem of the so-called Skorniakov Ter- 
Martirosian equation [27]] ■ For three identical bosons in 
2D where IZ 1 ^ 2 = TZ 1 ^ 3 = TZ 2 ± ^ 3 . a similar eigenequa- 
tion has been already derived in [U [12, (in that case, 
i] = 2 and y = 1/2 in Eq. (|2"2l ). An analogous integral 
equation has been obtained in the context of quasi-2D 
heteronuclear trimers in Ref. [29j where the third direc- 
tion is also taken into account. 



III. 1 + 2-BODY BOUND STATES 

In this last part, we show that Eq. (|2"2")l supports the 
existence of trimers which can be considered as shallow in 
the regime defined by Eq. (jTHJ) . We denote their binding 
energy by E tlim : 



E 



-Et 



< 



(27) 



The kernel in Eq. \22\ has the cylindrical symmetry, 
hence trimers states can be labeled by a radial quan- 
tum number and an orbital quantum number denoted 
below by m. We thus expand the source amplitude on 
the m-partial waves: 1Z(k) = J2m=o cos(m9)TZ m (k) , with 
9 = Z(e x ,k). The angular integration of the kernel can 
be performed for each m-partial wave, and one obtains 
the following integral equation: 



n m (k) r)(-iy 



/2D(«k) 



2ky 



00 It 



VT 2 



1 



Vt 2 



-n m (u), (28) 



where t = (u 2 + k 2 + q 2 )/(2yku) > 1. Eqs. (|24l27l) imply 
that /2D(«k) = 1/ hi(Ca 2 D|Kk|) is positive. Whence, for 
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i] = +1 (atoms of mass M are bosons) Eq. ((28)) admits 
solutions for to even only, while for r\ = — 1 (identical 
atoms are fermions) the orbital quantum number of the 
trimers is necessarily odd. We have computed the spec- 
trum for mass ratios x up to 200 -such high values can 
be reached experimentally using an optical lattice [14| . 

Bosonic spectrum- The bosonic spectrum is plotted 
in Fig. ([T]). For x < 200, we have found trimers with 
to < 8 [30| and the deepest trimer is in the m = sector. 
Apparition threshold of three first excited trimers in the 
m = sector are located at x = 1.77, x — 8.34 and x — 
18.27. Trimers of increasing orbital momentum appear 
as x increases: as an example, m = 2 trimers appear for 
x > 12.68. An important feature of the bosonic spectrum 
is that for all possible values of the mass ratio, there 
exists at least one trimer state. In the specific case where 
m i = M and to = 0, we recover the result of Ref. plj ]: 
there is only one bound state with a binding energy at 
E 3 /E dim = 2.36 • • • 




FIG. 1: Trimers' spectrum for two identical bosons of mass 
M interacting with a third atom of mass mj. The bind- 
ing energies Etrim are expressed in unit of the heteronuclear 
dimer's energy -Bdim in Eq. (|14[) . as a function of the mass 
ratio x — M/rm. One can notice that away from threshold, 
the bound state energy of a m-trimer is quasi-degenerate with 
the energy of a trimer with a lower m. 

Fermionic spectrum- The fermionic spectrum corre- 
sponds to the case where the two identical atoms are 
fermions (rj — — 1) and is plotted in Fig. ([2]). We find that 
the deepest shallow trimer has a p-wave internal angular 
momentum (m = 1) and exists only for x > x\ ~ 3.33. 
For increasing values of x > x\, other shallow trimer 
states appear in the to = 1 sector and for higher odd 
values of to. 

Numerical values of mass ratios for the threshold of 
apparition of 1 + 2-trimers in the to < 8 momentum 
sectors are gathered in Tabs. (IIIIip . 

Discussion- One important implication of the exis- 
tence of these universal trimers concerns losses resulting 
from two-body collisions of heteronuclear shallow dimers. 




FIG. 2: Same as Fig. |T} for two identical fermions of mass 
M interacting with a third atom of mass rrn. 





m = 


m = 1 


m = 2 


m = 3 


m = 4 


m = 5 


m — 6 


n = 




3.340 


12.68 


26.89 


45.95 


69.83 


98.84 


n = 1 


1.770 


10.41 


23.76 


41.98 


65.02 


92.97 


125.8 


n = 2 


8.341 


20.85 


38.17 


60.36 


87.49 


119.3 


155.9 


n — 3 


18.27 


34.59 


55.91 


82.08 


113.1 


149.0 


189.7 


n = 4 


31.62 


51.76 


77.04 


107.0 


142.09 


181.8 




n — 5 


48.33 


72.28 


101.5 


135.5 


174.2 






n — 6 


68.38 


96.23 


129.2 


167.0 








n = 7 


92.01 


123.5 


160.1 










n = 8 


119.0 


153.8 


194.8 










n = 9 


149.4 


188.2 












n — 6 


182.8 















TABLE I: Threshold values of the mass ratio M/rrii for the 
apparition of 1 + 2-body bound states. The momentum of 
the states is labeled by m while n is the hyperradial quantum 
number. Computations have been limited to mass ratios less 
than 200. 



In the case of bosonic identical atoms of mass M [see 
Fig. ([T])] the binding energies are -Etrim > 2i?dim for all 
values of the mass ratio. Thus, collisions of two shallow 
heteronuclear dimers can always lead to the formation 
of one {MMrrii )-trimer plus a single impurity atom of 
mass mi . Therefore, two-body losses prevent the exis- 
tence of a stable 2D gas composed of such dimers at equi- 
librium whatever the statistics of the impurity (atoms of 
mass TOi). For fermionic atoms of mass M, the deepest 
trimer has a binding energy -Etrim > SE^im for mass ra- 
tios a; > X2 — 18.3. Thus two-body collisions of heteronu- 
clear dimers where both atomic species are fermions, are 
elastic for a mass ratio x < X2 [31(- Whence a gas com- 
posed of such shallow dimers is expected to be relatively 
stable for x < x^. it is the case for 6 Li- 40 K or 6 Li- 87 Rb 
mixtures strongly confined in 2D providing Eq. (fl"6|) is 
satisfied, while for a 6 Li- m Y mixture the critical mass 
ratio is exceeded. In order to get some informations 
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m — 7 


m = 8 


n = 


132.3 


171.0 


n = 1 


163.5 




n = 2 


197.9 





TABLE II: Same as in Tab. for higher momentum. 
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FIG. 3: Plot of the probability in Eq. fl29j P(k < Q) at 
Q = 9dim as a function of the mass ratio x for the differ- 
ent trimers composed of bosonic atoms of mass M with the 
binding energies of Fig. {]]). Continuous lines: m = states, 
dotted lines: m — 2 states, dashed lines: m — 4 states. 



about this inelastic process, we consider for a given nor- 
malized trimer state the probability: 



P(k <Q)= f 



d 2 ki d 2 k 2 



(29) 



P(k < Q) in Eq. (|29|) gives the order of magnitude that 
the three particles of the trimer are separated each one 
from the others by distances greater than or of the order 
of 1/Q. The results for fermionic atoms of mass M and 
Q = <Zdim is given in Fig. Q and a similar plot for bosons 
is given in Fig. ©. These figures show clearly that the 
probability is maximum nearby the formation threshold 
of a trimer. Moreover, the probability at threshold de- 
creases rapidly as the value of the internal angular mo- 
mentum increases. In Fig. (J5J we have plot P(k < Q) as 
a function of l/Q for a m = 1 trimer near its threshold 
of apparition (x = 153.9): the probability drops out dra- 
matically as one considers increasing values of the inter- 
particle spacing: for Q = gdim/2, it is still greater than 
10% while for Q = Qdim/10, it is less than 1%. There- 
fore, we can conclude that the spatial extension of the 
lowest energy trimers is of the order of the dimer's size. 
Moreover, it has been shown that the long range effective 
heteronuclear pairwise interaction is repulsive for high 
mass ratio with a barrier located at ~ 1 /^dim [14| . Thus 
depending on the statistics, these results show that the 
dimer-dimer inelatic decays populate predominantly the 
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FIG. 4: Same as Fig. ([3]) for fermionic atoms of mass M with 
binding energies of Fig. (J2J. Continuous lines: m — 1 states, 
dotted lines: m — 3 states, dashed lines: m — 5 states. 
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FIG. 5: Plot of the probability in Eq. P(k < Q) as a func- 
tion of 1/Q in unit of 1/gdim for a m — 1-trimer at threshold 
(x = 153.9 and q - g d im < <7dim). 



to = or m = 1 shallow trimers nearby their apparition 
threshold (and also to less extents the to = 2 and m = 3 
trimers). Furthermore, from Figs. (|3l4j) it appears that 
such low energy 3-body bound states of size ~ 1/gdim ex- 
ist for all high mass ratios. Hence this motivate future 
detailed study of such four-body processes which are the 
most important source of 2D-losses even for high mass 
ratios. 



CONCLUSION 

In this paper, we have derived the low energy three- 
body wave equation using a ZRP model in the momen- 
tum representation for atoms strongly confined in a pla- 
nar wave guide. In this representation Green's functions 
are simply rational fractions, as a consequence techni- 
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cal implementation of the formalism is greatly enhanced 
with respect to the usual configurational space approach. 
Generalization to others dimensions and other few-body 
problems can be easily performed [24j . 

Depending on the statistics and on the mass ratio, we 
have predicted the existence of quasi-2D heteronuclear 
trimers composed of two identical atoms and of one an- 
other atom species. For a planar wave guide achieved 
by using a harmonic trapping along the tight direction 
(harmonic oscillator length a z ) these trimers can be ef- 
fectively considered as quasi-2D in a regime where the 
3D heteronuclear scattering length 0,30 is negative and 
I 3d| *C a z . While for fermions the heteronuclear inter- 
action is naturally negligible as a consequence of the Pauli 
principle, for identical bosons we have considered situa- 
tions where this channel of interaction is also negligible. 
Considering the shallow trimers as the most important 
source of decay in (heteronuclear-dimer)-(heteronuclear- 
dimer) scattering processes, we have shown that a dimer 
gas composed of fermionic atoms is stable with respect to 
binary collisions for a mass ratio M/m, < 18.3 while for 
bosons the gas is never stable. In the fermionic case, for a 
fixed value of the areal dimer density n compatible with a 
2D strong confinement (i.e. na 2 z <C 1), knowing that the 
2D heteronuclear atomic scattering length can be tuned 
to an arbitrarily large value for 030 < 0, then the 2D 
gas parameter na^ D can reach high value. This opens 
the possibility of achieving a stable 2D strongly corre- 
lated regime where one can expect also that for a mass 
ratio M/rrii > 3.33 the shallow trimer state(s) affect(s) 
the low energy gas properties. Practical achievement of 



such a low energy dimer gas in harmonic traps requires 
very high aspect ratios such that a z -C a2D *C On where 
aii is the harmonic oscillator length along the direction 
perpendicular to the tight confinement. 

Interestingly, similar results concerning shallow het- 
eronuclear trimers have been found in 3D [32j where 
the zero-range approximation is also performed and the 
Shrodinger equation is solved within the hyperspheri- 
cal formalism. In 3D, the critical mass ratio for the 
first p-wave shallow trimer is found at 8.18 (a larger 
mass ratio than in 2D). The universal shallow p-wave 
trimers exists also only for positive 3D scattering length 
(«3D > 0), while in 2D the analogous universal shallow 
trimer states exist for a2D 3> « z and thus a^D < 0. More- 
over in Ref. (33|, exact results have been found for ID 
universal heteronuclear trimers. A recent Letter [29j | has 
explored three-body physics in quasi-2D atomic planar 
waveguide. In this last reference the effective range of 
the 2-body interaction and also the transverse direction 
of the wave guide are both taken into account in the 
three-body problem. Eq. (1221) in our paper appears thus 
as a 2D limit of this approach. 
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